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ABSTRACT 

Observations of high-rcdshift supernovae (SNe) open a novel opportunity to study the 
massive star population in the early Universe. We study the detectability of superlumi- 
nous SNe with upcoming optical and near-infrared (NIR) surveys. Our calculations are 
based on the cosmic star formation history, the SN occurence rate, the characteristic 
colour and the light curve of the SNe that are all calibrated by available observa- 
tions. We show that 15-150 SNe up to z ^ 4 will be discovered by the proposed 
Subaru/Hyper Suprime-Cam deep survey: 30 deg^ survey with 24.5 AB mag depth 
in 2-band for 3 months. With its ultra-deep layer (3.5 deg^ with 25.6 AB mag depth 
in z-band for 4 months), the highest redshift can be extended to z 5. We further 
explore the detectability by upcoming NIR survey utilizing future satellites such as 
Euclid, WFIRST, and WISH. The wide-field NIR surveys are very efficient to detect 
high-redshift SNe. With a hypothetical deep NIR survey for 100 deg^ with 26 AB 
mag depth at 1-4 /im, at least ~ 50 SNe will be discovered at z > 3 in half a year. 
The number of the detected SNe can place a strong constraint on the stellar initial 
mass function or its slope especially at the high-mass end. Superluminous SNe at high 
redshifts can be distinguished from other types of SNe by the long time-scale of their 
light curves in the observer's frame, the optical colours redder than other corc-coUapse 
SNe and the NIR colours redder than any other types of SNe. 
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1 INTRODUCTION 

Core-collapse supernovae (SNe) are triggered by massive 
stars at the end of their lives. Thanks to the high luminosity, 
SNe can be detected in distant galaxies. By measuring the 
cosmic occurrence rate of core-collapse SNe, one can study 
the overall formation rate of the massive star population 
through cosmic time. With upcoming or planned large-scale 
optical and near-infrared (NIR) surveys, SNe could be de- 
tected to 2 ~ 2 or higher. It is thus important to derive a 
realistic estimate for the detectability in order to make an 
effic ient survey str ateg y for such observational programmes. 

iLien fc FieTd3 |2009l ) study the detectabil- 
ity of core-collapse SNe with the Large Synop- 
tic Survey Telescope (LSST , llvezic et al.l l2008l : 



ILSST Science Collaborations et al.l |2009| ). They show 



that the wide survey with LSST (with the 5 a r-band 
limiting magnitude of 25 mag) can detect Type II SNe, 
core-collapse SNe with hydrogen features, up to z ~ 2. 
Magnification by gravita tional lensing will bo ost the 
number of detectable SNe (|Oguri fc Marshallll2010l '). 



an extremely deep NIR survey ( 


Miralda-Escude & Reed 


ll997l:lMesineer, Johnson & Haiman 


l2006l~) with future satel- 
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lites such as James Webb Spac e Telescope (JWST). 
iMesinger. Johnson fc HaimanI (|2006l l show that up to sev- 
eral thousands Type II SNe can be detected at z; ~ 6 if 
the survey depth reaches the flux limit of 3 njy (30 mag in 
AB magnitude), which can be achievable with long (~ 10^ 
seconds) integration with JWST. 

It may be more promising to aim at detecting SNe that 
are bright in the rest frame optical to ultra-violet (UV) wave- 
lengths. Type Iln SNe, a type of SNe with narrow hydrogen 
emission lines, are ideal targets to detect at high-redshifts 
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l|Cooke|[200i ). Type Iln SNe are thought to be powered 
by the in teraction of th e SN ejecta with the circumstellar 
medium (ISchlege]|ll99d). a nd are bright in the rest frame 
optical to UV. ICooke et al.l |200i) reported the discovery of 
Type Iln SNe at z = 2.4 and 2.0 in the archival data of the 
legacy survey with Canada-France-Hawaii Telescope. 

There is another type of possibly luminous SNe 
worth mentioning here. Pair-instability SNe (PISNe) are 
thought to be triggered by very massive stars with 
mass in the ran^e of 140-260 Mq at zero metallicity 
(|Heger fc Woo slcv 200 ^), and are a lso expected to be very 



and SNe powere d by a magnetar (|Maeda et al.l l2007l : 
iKasen fc Bildsten.,20ia : IWoosIevl[201ol ) . 



luminous ( Scannapicco et al.ll2005l ). PISNe might also oc- 
cur even in the local U niverse, as suggested by the recent 
disco very of SN 2007bi (|Gal-Yam et al.ll2009l : iVoung et all 
I2OIOI ). However, PISNe are not very lumino us at rest 
UV wavelengths (|Kasen. Wooslev fc Hege"ill201lh . and thus 
they can be detected on ly at 2: < 2 even with LSST 
l|Pan. Loeb fc KasenI l201ll V It is possible, in principle, to 
detect PISNe at z > 6 with an extremely d eep NIR survey 
with K = 28.7 mag (|Scannapieco et al. |l2005l. see also recent 
papers by iPan. Kasen fc Loebl [20T1I : iHummel et"ai1l201lD . 



but the occurence rate of PISNe (or the formation rate of 
extremely massive stars) at such high redshifts is highly 
uncertain. Interestingly, recent theoretical studies and de- 
tailed simulations of the first generation of stars suggest 
that the characteristic mass of the first stars is not as large 
as the progenitor of a PISN (|Yoshida. Omukai fc Hernguistl 
I2OO7I : iHosokawa et al.l I2OIII ) . Then observing core-collapse 
SNe may be a more promising route to study the massive 
star population in the early Universe. 

Recently, s uperluminous SNe were discovered such 
as SN 2006gy lOfek ct al.' '2007"; 'Smith ct al.M 200l l2008l : 



lAgnoletto et al.l 12009 



Kawa bata ct al. 2009|) and SN 
(jQuimbv et al.l I2OO7I ) , where as other superlu - 
SNe include SN e 200 3ma (iRest et all l201lf) 
|2012| ). 2 008es jMiUer et al l 



2005ap 

minous 

2OO60Z (iLeloudas et al. 



I2OO9I : [Gczari ct all I2009D. 2008fz (|Drake et al.1 l2010t ). 
2008ma (Chatzopoulos et all \20n\). 201 0k x (PT FlOcwr 



IPastorello et al 



06F6 jBarbarv etal 
PTF09cwl 



2OIOI : iQuimbv et al.1 
' 2009j l 



l201ll ). SCP 

n— , PTF09atu, PTF09cnd, 

uimbv et al.1 I2OIID . PS l-lOky, PSl-lOaw h 



JChomiuk et al.ll201ll), and CSS10 0217 jPrake et al.ll201 j ) 



See also iRichardson et"aL ] (|2002D for other earlier candi- 
dates. These SNe have an absolute magnitude of about 
—22 mag in optical. Some of them are clas sified as Type 
Iln (e.g., lOfek et al.|[2007l : ISmith et aI.ll2007D . or Type IIL, 
which is a su bclass of Type II S Ne with the linear light 
curve decline (iMiller et al.1 12009| : ICezari et al] 120091). and 
others are hydrogen poor (Type Ibc, e.g.. IPastorello et al] 
l20ld : lQuimbv et al.|[201ll '). 

The nature of these superluminous SNe is poorly 
known. Some of them are powered by the circumstellar inter- 
action, or by the shock breakout from the dense circumstel- 
lar m edium (|Chevalier fc IrwinI I2OIII : iMoriva fc Tominagal 
I2OIII ) , as suggested by the presence of narrow emission lines 
in superluminous Type Iln SNe (e.g., [Agnolctto ct al..2009l ). 
It is also argued that superluminous SNe could be powered 
by a large amount of ^^Ni whic h is synthesized as a re- 
sult of energetic core-collapse SNe (lUmeda fc Nomotci|[200^ : 
lYoung et al.ll2010l : [Moriva et al.ll2010l '). Other scenarios in- 
clude the interaction b etween shells ejected by the pulsa - 
tional pair- instability (| Wooslev. Blinnikov fc Hegerl |2007^ ■ 



Although the origin of the high luminosity is not yet 
fully understood, the progenitors of the superluminous SNe 
are thought to be very massive. The long time-scale of the 
light curve is suggestive of a long diffusion times of opti- 
cal photons and hence of a large mass of ejecta (for a given 
kinetic energy). At least one progenitor of Type Iln SN is di- 
rectly identified to be as luminous as luminous blue variables 
with t he zero-age main sequence mass of Mzams > 5 — 
8OM0 (|Gal-Yam et allbOO?] : [Cal-Yam fc Leonardll2009^ . 



Superluminous SNe may provide an important infor- 
mation on the star formation in the high-redshift Universe. 
Because superluminous SNe are likely to be triggered by the 
very massive end of the stellar population, the relative oc- 
curence rate of such SNe with respect to the overall star 
formation rate can be a sensitive probe of the slope of the 
stellar initial mass function (IMF). Theoretical studies sug- 
gest that the shape of the IMF may be different in t he early 
Univ erse, or could even be of top-heavy one (e.g., iLarsonl 
Il998t ). Detecting superluminous SNe at high redshifts will 
give a direct evidence for such non-standard IMF. 

In this paper, we study the detectability of superlumi- 
nous SNe at high red s hifts with the upcoming wide-field sur- 
veys. iQuimbv et al.l (|201ll ) noted that superluminous SNe 
are bright enough that they can be detectable even at z ~ 4 
with 8m-class telescopes, such as Subaru and LSST. We de- 
rive a realistic estimate for the detectability for a few par- 
ticular set of surveys. For the optical survey, we focus on 
the pr oposed survey with t he Subaru/Hyper Suprime-Cam 
CHSC. lMivazaki et al.ll2006l '). The method of mock observa- 
tions is described in Section [2] We show that with a real- 
istic survey strategy we can detect superluminous SNe to 
2 ~ 4 — 5 (Section [3|. More emphasis is made on the de- 
tectability by upcoming NIR survey utilizing future satel- 
lites, such as Euclidj, The Wide-Field Infrared Survey Tele- 
scope (WFIRST) Q, and Wide-field Ima ging Surveyor for 
High-redshift (WISH) (Section li)). We study how the ex- 
pected number of detection at z > 3 is affected by the slope 
of the IMF at the massive end. Details of the sample selec- 
tion are described in Section [5| Finally, we give our conclud- 
ing remarks in Section [Sj 

Throughout the paper, we assume the Qm — 0.3, 
SIa =0.7 and Ho = 70 km s~^ Mpc~^ co smology, which is 
consistent with, e.g.. lKomatsu et al.l (|201ll ). The magnitudes 
are given in the AB magnitude when nothing is mentioned. 



2 SETUP FOR MOCK OBSERVATIONS 

We perform simulations of mock observations to study the 
detectability of superluminous SNe at high redshift. To this 
end, we need (1) spectral energy distributions (SEDs) of su- 
perluminous SNe and their time evolution, (2) cosmic occur- 
rence rate of superluminous SNe, and (3) a few key observa- 
tional parameters, i.e., the depth, the area, and the cadence 
of observations. In the following subsections, we describe the 
setup for our simulations. 



http: //sci ■ esa. lnt/6uclld| 
http : //wf irst . gsf c ■ nasa ■ gov) 



http : //www ■ wlshmlssion ■ org/en/lndex ■ html | 
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Rest wavelength (A) Days (rest frame) 

Figure 1. Rest- frame SED (left) and absolute light curves (right) of superluminous SN 2 008es, used as the b ase model of our simulations. 
The evolution of the luminosity and the effective blackbody temperature are taken from iMiller et al. 
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Figure 2. SEDs of hydrogen poor superluminous SNe PTFOQcnd 
(pink) and PTFlOcwr (SN 2010gx, green) compared with 
the blackbody wit h T = 15, 000 K. Spectral data are 
from i Q uimbv et alj 1I2OIII') and photomet ric data are from 
iQuimbv et alJ l l201ll) : |Pastorello et alJ l l2010h . 



2.1 SED of Superluminous Supernovae 

As our fiducial model, we use the obser ved SED of SN 2008es 
l|Miller et al.ll2009l : ICeiari et al.ll2009l ). SN 2008es is a well- 
observed object among superluminous SNe. The SED can be 
approximately described by blackbody radiation with tem- 
perature of 6000 — 15000 K depending on the phase. We 
adopt the model by iMiller et al.l (|2009l ) for the luminosity 
and the blackbody temperature. Note that the SED is not 
corrected for host galaxy extinction. We assume a cut-off of 
the flux at the wavelengths shorter than the Lyman limit. 
The model SED and rest-frame light curves for SN 2008es 
are shown in Figure [T] Around the maximum brightness, the 
blackbody temperature is ~ 15000 K, and the peak wave- 
length is found at about 3000 A. 

It is known that there is a variety of SEDs of super- 
luminous SNe; some SNe show H features while other do 
not. Figure [2] shows t he SEDs of hydrogen poor superlumi- 
nous SNe P TFOQcnd jQuinibv et al.ll201lh and PTFl Ocwr 
(SN 2010gx. |Pas"torello et al.ll2010l : lQuimbv et al.ll201ll ). For 
comparison, the blackbody spectrum with T = 15, 000 
K is also shown, which is assumed in our fiducial model 



at the maximum brightness. T he overall SEDs are simi- 
lar at the relevant wa velengths. IQuimbv" et al.1 (|201ll ') and 
IChomiuk et~all l|201lh also find that the evolution of the 
effective blackbody temperature is similar among the super- 
luminous SNe. Although the SEDs of SN 2008es may not 
apply to all the superluminous SNe, it is reasonable to as- 
sume that sup erluminous SNe hav e a bluer colour than other 
types of SNe l|Quimbv et al.ll201ll 'l. 

The characteristic time scale of the luminosity evolu- 
ti on differs among sup erluminous SNe (see e.g., Figure SI 
of IQuimbv et al.]|201ll ). For example, the bolometric light 
curve of SN 2008es declines by 0.85 dex in 50 days while that 
of PTF09cnd declines only 0.25 dex in the same period. Our 
fiducial model is based on SN 2008es, which yields conserva- 
tive estimates for the number of detections. We also perform 
simulations by using the blackbody model for PTFOQcnd by 
IQuimbv et al.l (|201ll ) to see the dependence on the model. 

We also note that dense observations of SN 2008es were 
done only after its maximum brightness. We do not extrap- 
olate the light curve before the maximum (Figure[T]), and do 
not count SNe in the early phase. Consequently, our simula- 
tions give conservative estimates. We argue that, if the light 
curve has a symmetric shape, the number of detections can 
increase up to by a factor of 2. 

SN 2008es has a peak magnitude of —22 mag in g-band, 
whereas some of superluminous SNe have even a brighte r 
magnitude up to -23 mag (e.g., IQuimbv et"aIll2007l . [2oTll '). 
It is not yet clear if the superluminous SNe are the objects 
at the luminous end of the luminosity function of normal 
core-collapse SNe, or they form an independent population. 
For our simulations, we assume the latter case, and adopt a 
Gaussian luminosity function with the g-band peak magni- 
tude of —22 mag. Since the dispersion of the peak magnitude 
is not well constrained, we adopt a small value (0.3 mag) 
so that the results are not affected by the luminous end of 
the luminosity function. Note that this luminosity function 
crudely includes the effect of extinction in the host galaxies, 
because the host extinction is not corrected in most of the 
literature. 

To perform realistic mock observations, we also include 
Type la, Type IIP, and Type Ibc SNe in our simulations. 
For the SEDs of these types, we use spectral templates by 
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iNugent.' Kim fc PerlmutteJ (|2002l ) Q. For their luminosity 
functions, we assume observed average peak i?-band Vega 
magnitudes and dispersion of nearby SNe, i.e., M_r = —18.6 
mag (cr =0.8 mag), —16.1 mag (cr =1.2 mag) and —16.1 
mag (c r =1.4 mag), for Type la, IIP and Ibc SNe, respec- 
tively (|Li et al.ll201ll ). 

2.2 Supernova Rate 

Superluminous SNe are likely to be the explosions of mas- 
sive progenitors. It is reasonable to assume that the occur- 
rence rate is proportional to the cosmic star formation rate 
p*{z) without significant time delay. We assume the mini- 
mum zero-age main sequence mass for superluminous SNe 
to be Mmin.sLSN = SO Mp) , which is the mini mum mass of lu- 
minous blue variables (| Gal- Yam et al ]|2007l V Suppose only a 
small fraction /slsn of massive stars (with > SOMq) explode 
as superluminous SNe. Then the rate of the superluminous 
SNe as a function of redshift -Rslsn(z) can be calculated as 



-RsLSn(2) = ./"sLSN p. (2) - 



J A. 



i){M)dM 



r 

J A, 



where tlj{M) is the stellar IMF 



(1) 



(2) 



We adopt a m o dified Salpeter A IMF of 
iBaldrv fc Glazebrookl (|2003l ') with the slope F = 0.5 
for O.lMo < M < O.5M0 and F = 1.35 for 
0.5Mq < M < IOOMq. The other parameters in Eq. 
lU are, the minimum mass of stars M^nin = 0.1A/q, the 
maximum mass of stars Mmax = IOOMq , and the maximum 
mass of superluminous SNe Mmax.SLSN = IOOMq. For the 
sta r formation rate, w e use t he parametrized prescription 
by iHopkins fc BeacomI (|2006h with a modified Salpeter A 
IMF. 

The fraction of superluminous SNe, /slsn, is not well 
determined. It is indeed one of the important goals of up- 
coming surveys to anchor the rate of such events to under- 
stand the population of progenitors. In the present paper, 

* |http : //supernova. Ibl . gov/- nugent/nugent_t6mplates .html | 



we cali brate this frac t ion w ith the available observations 
so far. iQuimbv et all (|201lf ) estimated the rate of hydro- 
gen poor superluminous SNe to be about 10~* Mpc~^ yr~^ 
at 2 ~ 0.3, which is only 10~* of total core-collapse SNe. 
This fraction corresponds to /slsn = 2 x 10~^ of massive 
stars with > 50Mq . Including hydrogen rich superluminous 
SNe, such as superluminous Type Iln SNe, the rate can be 
higher. Thus, in our simulations, we adopt /slsn = 2 x 10""^ 
- 2 X 10"^, i.e., 10"* - 10"^ of total core-collapse SNe are 
superluminous SNe. 

The left panel of Figure [3] shows the rate of the su- 
perluminous SNe. The rate lies roughly between 10~^ and 
10~^ Mpc~^ yr~^ in the redshift range we are most inter- 
ested in. The SN rate per volume per year can be readily 
converted to the number per redshift per field (solid angle 
Q) and per a certain duration of observations as 



diV 



dn dt dz 



1 dKom 

1 + z dQdz ' 



(3) 



where Vcom is the comoving volume. The right panel of Fig- 
ure [3] shows the number of SNe per 1 square degree, per 1 
year, and per redshift. The number of superluminous SNe 
at 2: ~ 2 is about unity within 1 deg^ per year, and does not 
decrease significantly to z ~ 4. 



2.3 Survey Strategy 

We primarily consider three types of surveys with different 
survey area, depth, cadence, and wavelength (Table [ij . In 
practice, deep imaging observations must be performed be- 
fore each survey to make a reference image. The transient 
survey described below will be performed after such refer- 
ence observations. 

First, we adopt the proposed parameters for the deep 
layer of the Subaru/HSC survey (hereafter called HSC- 
Deep). The survey is planned to cover 30 deg^ The same 
field will be visited twice per month with a 6-day cadence 
around the new-moon phase. Such observations will be re- 
peated for continuous 3 months. 

We define the survey cadence with 3 parameters; ni 
and At are the number of visits and its cadence within one 
month, and n2 is the number of monthly visits. Thus, for 
HSC-Deep, ni = 2, At = 6 days, and 712 = 3. 
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Table 1. Survey Parameters 



Survey 


Area 


At 


"1 


712 




5 (T limitinf 


5 magnitude per visit 






(dcg2) 


(day) 








rtlr 


mi 


m2 


my 


Subaru/HSC Deep 


30 


6 


2 


3 


26.1 


25.8 


25.6 


24.5 


23.2 


Subaru/HSC Ultra Deep 


3.5 


6 


3 


4 


26.9 


26.6 


26.6 


25.6 


24.3 














'"F200W 








NIR Deep 


100 




1 


6 


26.0 


26.0 


26.0 


26.0 


26.0 



The HSC survey is proposed to be performed with 5 
broad band filters [g,r,i,z, and y). The planned survey 
depth is summarized in Table [T] For HSC-Deep, the expo- 
sure time required to reach the designated depths is 12 min 
for g and r bands, and 18 min for i,z, and y bands. 

The second survey we consider is the Ultra Deep layer of 
Subaru/HSC survey (hereafter HSC-UltraDeep). It will ob- 
serve two Subaru/HSC fields (3.5 deg^ in total) with deeper 
magnitude limits and more frequent cadence (see Table [T|. 
The required exposure times for these limits are 30 min, 30 
min, 60 min, 90 min, and 90 min for g, r, i, z, and y band, re- 
spectively. The same field will be observed 3 times per month 
with a 6-day cadence, and this sequence will be repeated for 
4 months, i.e., ni = 3, At = 6 days, and n2 = 4. 

Finally, we explore SN detection by NIR surveys. We 
propose a hypothetical survey with a wide area of 100 deg^, 
which will be made possible by future satellite missions such 
as Euclid, WFIRST and WISH. For the transmission curve 
of NIR filters, the broad band filters of JWST are used 0. 
We assume a constant depth of 26.0 mag for 5 bands in 1-5 
/^m. For the NIR survey, we assume only one visit per month 
(ni = 1) for continuous 6 months (n2 = 6). 



2.4 Methods of Simulations 

In this subsection we describe the details of our mock obser- 
vations. We first generate superluminous SNe according to 
the SN rate in the right panel of Figure [3] and the adopted 
survey area assumed for each survey strategy. The observed 
light curve of each superluminous SN at a certain redshift 
is calculate d from the SEP in Figure [T] with appropriate K- 
correction (|Hogg et al.ll2002l ). Here the luminosity of each 
SN is assigned according to the luminosity function (with a 
Gaussian probability distribution). The date of explosion is 
determined by a random number. In this way, we generate 
expected hght curves of superluminous SNe at a range of 
redshifts. 

Then we perform mock observations. The magnitude 
of superluminous SNe are compared with the limiting mag- 
nitude of the survey specified by At, ni and n2 (see Table 
[1]). We do not consider contamination of the host galaxy, 
i.e., the detection limit in Table [T] is kept constant for all 
the objects. Strictly speaking, this is not the case when an 
object is discovered near a bright galaxy. However, for high- 
redshift superluminous SNe, the host galaxies are expected 



^ |http: //www. stscl ■ edu/jwst/instruments/nlrcam/] 
|lnstruinentdesign/f llters/lndex_htm'T] 
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Figure 4. Detection efficiency of superluminous SNc as a function 
of redshift for HSC-Decp (green) and HSC-UltraDeep (blue), and 
for a hypothetical NIR deep survey (red). 



to be faint (Section [5|, and thus it is a sound approximation 
to neglect the contamination. 

We impose stringent detection criteria. Only SNe that 
fulfill both of the following two criteria are counted as "de- 
tected"; (1) brighter than 5 a detection limit in more than 
2 bands at least at one epoch, and (2) brighter than the 5 cr 
detection limit at more than 3 epochs at least in one band. 
Typically, we generate hght curves for 1000 superluminous 
SNe per one run. The simulations are performed 1000 times, 
and the number of the detection is averaged over the real- 
izations. 



3 DETECTABILITY WITH OPTICAL SURVEY 

We first calculate the detection efficiency for the three sur- 
veys we propose. The detection efficiency is defined to be 
the ratio of the detected number of SNe to the total num- 
ber of SNe exploded in the survey period. The green and 
blue lines in Figure |4] show the detection efficiency (/detect) 
of superluminous SNe for HSC-Deep and HSC-UltraDeep, 
respectively. We see rapid drop in /detect at z 3 — 4 and 
z ^ 4 — 5 for Deep and for Ultra Deep, respectively. These 
limiting redshifts corresponds to the redshifts where super- 
luminous SNe with the mean absolute magnitude —22 mag 
becomes undetectable in 2-band (see Appendix |X] for the 
SEDs in observer's frame). HSC-Deep is almost complete 
up to z ~ 3 while HSC-UltraDeep is complete up to z 4. 
The efficiency at lower redshift is not unity because of our 
somewhat stringent detection criteria. If we impose a looser 
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Figure 5. An example of the simulated t-band liglit curve of a 
superluminous SN (red) for HSC-Deep. The light curve is com- 
pared with that of Type la (black), Type IIP (blue), and Type Ibc 
(green) SNe at lower redshifts. Photometric errors shown here are 
computed according to the limiting magnitudes in Table [T] The 
light curves of Type la SNe and Type Ibc SNe are shifted toward 
the right for clarity. 



criteria, for example, > 5(t detection at only 2 epochs, the 
efSciency becomes nearly unity at low redshifts. The effi- 
ciency and its redshift dependence are not largely affected 
even when PTF09cnd, which has a slower hght curve than 
our fiducial model, is used as our input model. 

It is worth mentioning that SNe exploded before the 
survey period can also be detected. This effect is not in- 
cluded in the detection efficiency, but such SNe are included 
in the detection number shown below. 

Figure [5] shows a few examples of simulated light curves. 
Photometric errors are computed according to the limiting 
magnitudes in Table [T] The i-band hght curve of a super- 
luminous SNe is shown by the red points. It is compared 
with those of other types of SNe. Superluminous SNe at 
2-^3 have comparable magnitude with other types of SNe 
at 2 < I. Since the observed i-band wavelength corresponds 
to the rest-UV wavelengths, the decline rate of superlumi- 
nous SNe is not extremely long. The observed time-scale is 
slightly longer than that of Type la and Ibc SNe, but shorter 
than that of Type IIP SNe at plateau phase. 

When PTF09cnd is used as input, the observed bright- 
ness of a model superluminous SN at z = 3 decreases as 
slowly as that of Type IIP SNe sd, z = 0.3. In this case, about 
50 % of the superluminous SNe detected with the HSC-Deep 
do not show the variability larger than 1 mag in any band. 
It is thus important to have deep reference images before 
the survey, otherwise such SNe would be missed. 

The expected number of detection of superluminous 
SNe with each layer of the survey is shown in Figure [HI 
Two solid lines show the cases for /slsn = 2 x fO^^ (top) 
and 2 x 10"'^ (bottom). Our reasonable guess is that the 
true number will lie in between them. The total number of 



^ The fraction of SNe with a small variability (< 1 mag) is only ~ 
15 % when the model of SN 2008cs is used. This fraction becomes 
smaller if the baseline the survey period is longer. For example, 
if the HSC-Deep is performed over 6 months keeping the total 
number of visit, i.e., ni = 1 and n2 = 6, the fraction of the small 
variability case is only 20 % for the PTFOQcnd model, and almost 
% for the SN 2008es model. 
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Figure 6. Expected number of detection as a function of redshift 
for HSC-Dccp (upper panel) and HSC-UltraDeep (lower panel). 
The two solid lines show the cases for /slsn = 2 x 10^'^ (top) 
and 2 X 10^"^ (bottom). The total number of detection is 15-150 
for HSC-Deep and 3-30 for HSC-UltraDeep. 



detection is 15-150 for the Deep layer, and 3-30 for the Ul- 
tra Deep layer. If PTF09cnd is used as input, these numbers 
increase by 20 % because of the slow luminosity evolution. 

For both the layers, the median redshift is found to be 
z ~ 2, where the star formation rate and hence the intrin- 
sic SN rate has a peak. The wider survey area of HSC-Deep 
gives a larger number of detection sd, z — 1 — 3. In the future, 
a dedicated deep survey by LSST will discover more super- 
luminous SNe. If a simi lar depth is achieved, for example, 
in th e deep drilling field (|LSST Science Collaborations et al.l 
l2009t ). the expected number of detection is simply propor- 
tional to the survey area. 

We expect that the rate of superluminous SNe at 
2: ~ 1 — 3 is accurately determined with HSC-Deep. HSC- 
UltraDeep can possibly detect superluminous SNe even at 
2~5if/sLSN = 2xlO-^This is in contr ast to PISNe, which 
we ca n detect only at z < 2 with LSST ((Pan. Loeb fc KasenI 
l201lf ). The highest redshift of superluminous SNe depends 
on the intrinsic rate of superluminous SNe, which can be ob- 
servationally constrained by HSC-Deep. We emphasize that 
surveys with the 2 layers are important to understand first 
the properties of the superluminous SNe (Deep layer) and 
then to detect highest-redshift SNe (Ultra Deep layer). 
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Figure 7. NIR light curves of superluminous SNe (red) com- 
pared with Type la (black), IIP (blue), and Ibc (green) SNe. 
For superluminous SNe, Type la and IIP SNe, magnitudes in the 
F200W filter is shown. Since the SED of Type Ibc SNe at 2^m 
is uncertain (see Figure fAll l. magnitudes in the F115W filter is 
shown for Type Ibc SNe. The light curves of Type la SNe and 
Type IIP SNe are shifted toward the right for clarity. 
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Figure 8. Expected number of detection as a function of redshift 
for the 100 deg-^ NIR survey. The two solid lines show the cases 
for /sLSN = 2 X 10-2 (top) and 2 X IQ-^ (bottom). The total 
number of detection is 120-1200. 




Figure 9. Expected cumulative number of superluminous SNe at 
2 > 3 as a function of the sloop index P of the IMF at M > 50Mq . 



For this plot, our 'conservative' 
assumed. 



SN rate (/slsn = 2 X 10"^) 



4 NIR SURVEY 

4.1 Expected Number of Detection 

It is likely that a NIR survey is more efficient than the optical 
surveys to detect high redshift superluminous SNe. The red 
line in Figure [4] shows the detection efficiency for our pro- 
posed NIR survey. The efficiency is approximately constant 
at z = — 6. The maximum efficiency is 0.5 simply because 
we require detection at 3 epochs. Most superluminous SNe 
up to z ~ 6 are above the detection limit. However, because 
of our criteria, superluminous SNe occured in the latter 3 
months in the survey period are not included in Figure S) 
Using NIR bands is essential to perform a complete survey 
at redshift z > 3. 

Figure [7| shows the simulated light curves of superlumi- 
nous SNe (red) compared with other types of SNe. Thanks 
to the high rest-UV luminosity (Figure [T]), superluminous 
SNe even at z = 6 can be as bright as 24-25 mag at 2 fim. 
Because of the time dilation, the fight curve of superlumi- 
nous SNe at high redshifts has a long time-scale in NIR. 

The expected number of detection is shown in Figure [8] 
It is highly interesting that, with a 100 deg^ survey, hundreds 
of z > 2 SNe can possibly be discovered. The total number 
of SNe is 120 and 1200 for the cases of /slsn = 2 x 10"^ 
and 2 x 10"^, respectively. We discuss the usefulness of such 
a large sample at high redshifts in the next subsection. 



4.2 Sensitivity to the IMF 

As we mentioned earlier, the superluminous SNe of our tar- 
get are thought to be triggered by very massive stars. The 
occurrence rate is very sensitive to the number fraction of 
such massive stars. This in turn can be used, if the occurence 
rate is accurately d etermined, to prob e the slope of the IMF 
at the massive end (jCooke et al.ll20()9l ) . We calculate the sen- 
sitivity of the expected number of the detection to the slope 
of the stellar IMF. For simplicity, we fix the minimum mass 
of superluminous SNe (Afmin.sLSN) and the fraction of su- 
perluminous SNe (/slsn) for the following comparison. Note 
that these parameters can be observationally constrained by 
the optical survey as we discussed in Section 3. To study the 
influence of the relative fraction of massive stars, we change 
the slope of the IMF at the large mass > 50Mq, whereas 
keeping the slope at < 5OM0 to be 1.35. It is important to 
note that changing the IMF at the high-mass end has lit- 
tle impact to the calibration of the star formation rate since 
most UV continuum emission comes from massive stars with 
< 50Mq. Thus the occurence rate of superluminous SNe is 
a robust probe of the stellar IMF. 

Figure |5] shows the cumulative number of detection at 
z = 3 — 6 as a function of F, the slope of the IMF. Here the 
fraction of superluminous SNe is conservatively assumed to 
be /slsn = 2 X 10"''. Even in this case, with F < 1.1, the 
expected number can be about three times as much as that 
with F = 1.35. Thus, if the star formation rate is indepen- 
dently measured, and if the SN rate can be determined with 
an accuracy within a factor of 3, the very flat slope of the 
massive end of the IMF can be inferred, or rejected. 
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Figure 10. Optical (left) and NIR (right) colour-colour diagram for various types of SNc. Our template superluminous SN 2008cs at 
high rcdshifts (red) are compared with Typo la SNe (black), Type IIP SNe (blue), and Type Ibc SNe (green) at low redshifts. The 
redshift of Type la, IIP, and Ibc SNe arc selected so that the observed magnitudes are similar to those of high-redshift superluminous 
SNe. Colours at different epochs are connected with lines, and the square shows the colour at maximum brightness. In optical, colours 
of superluminous SNe at high redshifts are red both in r — j and i — z, compared with Type IIP SNc and Type Ibc SNe. Only Type la 
SNe at redshift 2 = 1.5 can have colours redder than superluminous SNe both in r — i and i — z. In NIR, superluminous SNe at high 
redshifts are redder than any other types of SNe with similar observed brightness. For the detailed SED, see Appendix 1X1 



5 NOTES ON THE SAMPLE SELECTION 

So far we have discussed the number of detectable SNe of a 
particular type, but have not fully evaluated the efficiency of 
target selection for actual surveys. Because superluminous 
SNe are intrinsically rare events, upcoming optical surveys 
using 8m-class telescope, such as Subaru and LSST, and 
NIR surveys with Euclid, WFIRST, and WISH, will discover 
more Type la and ordinary core-collapse SNe. These normal 
SNe will be more abundant than superluminous SNe by a 
factor of about 100. Thus, to identify superluminous SNe at 
high redshifts, we must select the candidates properly and 
effic iently. 

ICookd I 20081 ) propose a method to select Type Iln SNe 
from colour-selected Lyman break galaxies. The advantage 
of the method is that the photometric redshift of the host 
galaxy is already known. When Lyman break galaxies at 
z > 2 are selected, most of SN candidates found in this 
sample are likely to be Type Iln SNe or superluminous SNe. 

Here we consider the characteristics of superluminous 
SNe in a more general case. An important issue would be 
that the host galaxy of a distant SNe is likely to be too faint 
to be detected. We summarize the properties of (1) hght 
curve and its time scale, (2) colours, and (3) host galaxies 
of superluminous SNe at high redshifts. 

The observed time-scale is determined by the intrin- 
sic time-scale and the time dilation due to cosmological 
redshift. Superluminous SNe observed so far indeed tend 
to show an intrinsically long time-scale. For an extreme 
case, the brightness of SN 2006gy stayed within 1 magni- 
tude below the max imum for about 100 days. Other ob - 
jects such as 2008es dMiller et al II2009I : iGezari et al.ll2009h . 
SCP 06F6 (iBarbarv et a,l.l |2009|) . PTF09atu, PTF09cnd, 
PTFOOcwl iQuimbv et aLlbOllh also have a long time-scale. 
Their decline tends to be faster in bluer bands (Figure [I}. 

In optical surveys, we observe high-redshift superlumi- 
nous SNe at the rest-UV wavelengths. The observed time- 



scale of superluminous SNe is slightly longer than that of 
Type la and Ibc SNe, and shorter than Type IIP SNe (Fig- 
ure [SJ, when we use SN 2008es as a template. When we use 
PTF09cnd as a template instead, the light curve is as slow 
as Type IIP SNe. With such a slowly evolving light curve, 
it may be difficult to distinguish superluminous SNe from 
Type IIP SNe. 

In NIR surveys, the long time-scale of superluminous 
SNe is pronounced (Figure[7]). Since NIR wavelengths traces 
the rest-frame optical at 2: > 2, superluminous SNe at 2 > 2 
can stay detectable for more than 100 days thanks to both 
the intrinsic brightness and the dilation effect (Figure [7]). 
As a result, the time-scale is longer than any other types of 
SNe. 

The second diagnostic is the observed colour. Figure 
[TOl shows optical r ~ i vs. i - z (left) and NIR F115W- 
F200W ([1.2]-[2.0]) vs. F200W-F277W ([2.0]-[2.8], right) 
colour-colour diagrams. Colours of our template superlu- 
minous SN 2008es at high redshifts are shown in red. For 
comparison, we also show the expected colours of Type la 
SNe at 2 = 1.0 and 1.5 (black). Type IIP SNe at (blue), 
and Type Ibc SNe (green) aX z = 0.3 and 0.5. These SNe 
are expected to have similar observed magnitudes to those 
of superluminous SNe at z > 3. The colours of these types 
of SNe are calculated by using t he spectral templates by 
iNugent. Kim fc PerlmutteJ (|2002h . See Appendix |X] for the 
SEDs. The reddening in the host galaxy is not taken into 
account here. It is worth mentioning that there is a scatter 
in optical colour of observed SNe at low redshifts. Including 
the contribution from the host galaxy, the scatter is typically 
about 1.0 mag for Type IIP and 0.5 mag for Ty pe Ibc SNe 
(e.g., lOlivares E. et al.ll2010l : iDrout et al.ll201ll ). For Type 
la SNe, a wide color spread is known, bu t most objects hav e 
similar colours within 0.2-0.4 mag (e.g., IWang et al.ll200^ ). 

Superluminous SNe at high redshifts tend to have red- 
der colour than Type IIP and Ibc SNe both in r — i and i~z. 
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For superluminous SNe at high redshifts, the observed r,i,z 
wavelengths correspond to the near-UV wavelengths in the 
rest frame, which are at the blueside of the spectral peak. 
On the other hand, most of Type IIP and Ibc SNe with opti- 
cal surveys will be discovered at 2 < 1. Their observed r — i 
and i — z colours still trace the optical wavelengths in the 
rest frame. Thus, colours of superluminous SNe at high red- 
shifts tend to be redder than those of Type IIP and Ibc SNe 
(see also Appendix [A} ■ In optical wavelengths, only Type 
la SNe at z > 1 have a similarly red or redder colour than 
superluminous SNe.Q 

The red colour of superluminous SNe is more prominent 
in NIR wavelengths. See the right panel of Figure 1101 NIR 
wavelengths traces the blue optical or near UV wavelength of 
superluminous SNe at high redshifts (see Appendix |A]|. For 
other types of SNe, observed NIR wavelengths correspond 
to the red optical (Type la) or NIR (Type IIP and Ibc) in 
their rest frames. Therefore, NIR colours of superluminous 
SNe at high redshifts tend to be redder than those of other 
types of SNe. The red colour of superluminous SNe results 
simply from the fact that they are at high redshifts. Thus, 
their red color does not depend largely on the choice of the 
template SEDs. 

The last characteristic is the faintness of the host galax- 
ies of superluminous SNe. An obvious reason for this is 
the di stance to the host galaxies. Interestingly, iNeill et"al] 
l|201lf ) argue that the host galaxies of superluminous SNe 
have a faint intrinsic luminosity; most of them have an ab- 
solute r-band magnitudes of > — 20 mag[f|. They also have 
quite blue UV-optical colours NUV — r = 0. If there is no 
redshift evolution in the host galaxy population, the opti- 
cal magnitudes of the host galaxies at z > 3 are likely to 
be close to the detection limit. This is in contrast to Type 
la SNe at z = 1 — 1.5, whose host galaxies w ere identified 
in op tical wavelengths with m; = 23-26 mag (|Graur et al.l 
|20l3). Overall, faintness or non-detection of the host galax- 
ies of the high redshift superluminous SNe may indeed be of 
help with sample selection in addition to the characteristic 
features in the light curve and colour. 



6 CONCLUSIONS 

We study the detectability of superluminous SNe with up- 
coming optical and NIR surveys. By assuming that the pro- 
genitors of superluminous SNe are more massive than 50 
Mq and a fraction /slsn = 2 x 10"^ - 2 x 10"^ of these 
massive stars explode as superluminous SNe (or the rate of 
superluminous SNe are 10"* - 10"^ of total core-collapse 
SNe), the cosmic occurrence rate of superluminous SNe is 



about 0.5-5 deg ^ yr ^ redshift ^ at z ~ 1 — 3 and about 
0.1-1 deg"^ yr"^ redshift"^ at 2 ~ 4 - 5. 

We predict that 15-150 superluminous SNe up to 2 ~ 4 
will be detected with the Deep layer of the upcoming Sub- 
aru/HSC survey if all the objects meeting our detection cri- 
teria are indeed selected. About the half of them are at 
2 > 2. By the Deep layer, the cosmic occurrence rate of 
superluminous SNe at 2 ~ 2 — 3 can be anchored. The Ultra 
Deep layer will discover a smaller number of SNe (3-30), but 
SNe at 2 ~ 5 could be detectable. A similarly deep, ded- 
icated survey by LSST will discover more superluminous 
SNe, whose number is simply proportional to the survey 
area. 

We show that future NIR surveys will detect more than 
a few hundreds of superluminous SNe in half a year if 100 
deg^ field is surveyed with the depth of 26 mag. Even with 
the conservative estimate of the rate (2 x 10"'^ of massive 
stars with > 5OM0, or 10^** of total core-collapse SNe), 
such NIR survey will detect about 50 superluminous SNe at 
2 > 3. 

Superluminous SNe at high redshifts can be selected 
in the sampl e of high-redshift galaxies as demonstrated by 
ICookd (|2008l '). Also in the general survey, covering all red- 
shift ranges, superluminous SNe can be distinguished by 
other types of SNe by (1) the observed long time-scale espe- 
cially, (2) the optical colour redder than Type IIP and Ibc 
SNe, or the NIR colour redder than any other types of SNe, 
and (3) faint host galaxies. 

Detection of high-redshift SNe opens an exciting oppor- 
tunity to study the massive star population at high redshifts. 
The progenitor of superluminous SNe is likely to be very 
massive. Therefore, the detected number is very sensitive to 
the IMF. If the slope of the IMF is F < 1.1 at the high 
mass end (M > 5OM0), the detected number will be in- 
creased by a factor of >3. If the star formation rate at 2 > 3 
is measured with a small uncertainty by other methods, the 
detected number of superluminous SNe with the NIR survey 
can be an indicator of top-heavy IMF at the high-redshift 
Universe. 
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^ It is worth mentioning here the expected optical colours of 
quasars as possible contaminations. Most quasa rs at z < A have 
r — i = —0.5 - 0.5 and i — z = —0.5 - 0.5 l lSchneider et al] 
|2007|) . which are bluer than superluminous SNe at high redshifts. 
Quasars at 2 = 4 — 5 have a red r — i colour (up to ~ 2), which 
is comparable to the colour of superluminous SNe and Type la 
SNe. 

* Note however that there could be an observational bias that 
superluminous SNe are easily detected if their host galaxies are 
faint. 
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F277W, F356W, and F444W for NIR). In this figure, we as- 
sume peak y-band Vega m agnitudes that are brighter than 
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Figure Al. SEDs of various types of SNc at maximum brightness at different redsiiifts. SEDs of superluminous SNc (red) are compared 
witli other types of SNe at most probable redshift ranges. 
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